Bouncing ball tunneling in quantum dots 
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We show that tunneling through quantum dots can 
be completely dominated by states quantized on stable 
bouncing-ball orbits. The fingerprints of bouncing-ball tun- 
neling are sequences of Coulomb blockade peaks with strongly 
correlated peak height and asymmetric peak line shape. Our 
results are in agreement with the striking correlations of peak 
height and transmission phase found in recent interference 
experiments with quantum dots. 
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Transport experiments with semiconductor quantum 
dots yield important information on the energy levels 
and the wave functions of confined many-electron sys- 
tems |IJ. In the Coulomb blockade regime current can 
flow only if two different charge states of the quantum 
dot are tuned to have the same energy. This produces 
large nearly equally spaced conductance peaks as a func- 
tion of the dot potential (tuned via electrostatic gates) . 
At low temperatures, the magnitude of each conductance 
peak directly measures the magnitude of a single reso- 
nant wave function near the contacts to the leads. Some 
years ago, a statistical theory for the peak height distri- 
bution was developed based on the assumption that 
the wave functions are completely uncorrelated and de- 
scribed by random matrix theory. Recent experiments 
p,0 found good agreement with the predicted distribu- 
tions. At the same time, one of these experiments re- 
ported large correlation between the heights of adjacent 
peaks. Much stronger correlations of both the peak height 
and the phase of the transmission amplitude through a 
quantum dot were observed in novel interference experi- 
ments |^||. The correlations in these experiments com- 
prised all Coulomb blockade peaks found in a gate volt- 
age scan, yielding sequences of more than 10 correlated 
peaks. 

Theoretical work |Q demonstrated that nonzero tem- 
perature can partially account for but not fully explain 
the peak correlations observed in Ref . M . More recently 
PL it was shown that additional correlations can arise 
from short periodic orbits. Such orbits lead to deviations 
from random matrix predictions even for dots whose clas- 
sical dynamics is completely chaotic. The predicted peak 
modulations M are of the type observed in Ref. B , but 
can not account for the much stronger correlations in the 
interference experiments HH- Despite of considerable 
theoretical efforts f§-[ll]] the origin of these correlations 
has not been understood. In particular, the most strik- 



ing feature of the experiments , that the correlations 
comprise all observed Coulomb peaks, has not been ex- 
plained. We address this feature below. 

We show that tunneling through quantum dots can 
be completely dominated by wave functions quantized 
on stable bouncing-ball orbits. Bouncing-ball tunneling 
(BBT) dominates the transport if (i) the leads are at- 
tached opposite to each other and (ii) if they are suf- 
ficiently wide providing near normal injection of elec- 
trons into the dot. Under the conditions (i) and (ii) the 
dot conductance exhibits sequences of strongly correlated 
Coulomb peaks. All peaks within a sequence are domi- 
nated by the same bouncing-ball mode. The unique fin- 
gerprint of BBT is a characteristic asymmetry of the peak 
line shape for the peaks in the tails of the sequences. This 
asymmetry is caused by the breaking of the particle-hole 
symmetry as the bouncing-ball state moves away from 
the Fermi energy. We find the asymmetry in the line 
shapes measured in Ref. || , and suggest that the inspec- 
tion of line shapes be used as a diagnostic for BBT in 
other experiments on quantum dots. 

Our starting point is the transmission probability T r i ng 
through an Aharonov-Bohm ring with a quantum dot 
embedded in one arm. We assume that the dot is in the 
Coulomb blockade regime. The tunnel barriers around 
the quantum dot suppress multiple reflections of electrons 
across the ring, and T r i ng is given by pd| ] 

df 
dE, 



dE 



|to + i(£)exp[27ri$/$o]r. (!) 



Here, $ is the magnetic flux through the ring, $o is the 
flux quantum, and to, t(E) are the transmission ampli- 
tudes through the free arm and the arm with the quan- 
tum dot, respectively (to is assumed to be independent 
of injection energy). The ring is connected to two reser- 
voirs labeled up (U) and down (D) occupied according 
to the Fermi distribution /. Below, we investigate the 
phase (j>QD = arg J dE(—df /dE)t(E) and the transmis- 
sion coefficient T QD = / dE(-df/dE)\t(E)\ 2 . Both can 
be measured in quantum dot interference experiments 
ii- 

The transmission amplitude t(E) is related to the re- 
tarded Green function G pq of the quantum dot, 

t(E)=J2V p U V q D G pq (E), 

pq 



(2) 



where V p is the matrix element for tunneling between 
level p in the dot and the reservoir I = U,D, For a 



1 



weakly coupled dot the width of the level p is given by 



1 p 
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where T 



V; 



for decay into the reservoir I. 
expressed M as 



p \ 2 is the partial width 
The matrix element is 



V p = l^J / dsMs,z)d z r p (s,z) 



(3) 



where the integration is performed along the edge be- 
tween the potential barrier and the quantum dot (d z de- 
notes the derivative normal to the barrier). The wave 
function ip p corresponds to Dirichlet boundary conditions 
in the dot, while the barrier tunneling is fully included 
in the lead wave function tpi. The transverse potential in 
the tunneling region can be taken quadratic Q yielding 
ipi ~ c/exp[— (s — Sz) 2 /2cig ff ], where s is the transverse 
coordinate, si the center of the constriction and a e g its 
effective width. We can restrict ourselves to the low- 
est transverse mode since higher modes are suppressed 
by the barrier penetration factor (included in q). The 
Gaussian form of the lead wave function is convenient 
but not crucial for the results presented below. 

We obtain the dot wave functions tp p for an effec- 
tive potential which accounts for the dot confining po- 
tential and the interactions in the dot in a mean-field 
sense. We use a billiard approximation and model the 
dot by a hard-wall potential at the boundary, parame- 
terized in polar coordinates by R(4>) = Rq[1 + ecos(2</>)]. 
The parameter e measures the quadrupolar deformation 
out of circular shape. We use nonzero e to mimic |l2] ] 
the dots used in As in the experiments, we as- 

sume that the leads are attached opposite to each other 
at the boundary points closest to the origin (the points 
with (j> — ±7r/2). In contrast to most studies of quan- 
tum dots which assume the dynamics in the dot to be 
either regular or fully chaotic, our boundary parameter- 
ization yields mixed classical motion. This is illustrated 
for e = 0.2 in Fig. 2, showing chaotic motion in most of 
phase space except for two large islands associated with 
stable bouncing-ball motion between the contacts to the 
leads. We note that billiards of quadrupolar or similar 
shape have recently been used in studies of optical mi- 
croresonators [p^| , ^| . 

In Fig. 1 we show the coupling strength g p = 
T%Tg/(T%+Tg) over a sequence of 100 billiard states for 
two different values of a e g. Note that g p is proportional 
to the conductance peak height G p = (e 2 /h)(ir/2kT)g p 
measured in low temperature Coulomb blockade exper- 
iments [0. We calculated g p from Eq. (||) using quan- 
tum dot wave functions obtained numerically with the 
boundary integral method For narrow leads the 

peak heights fluctuate with the level index p, occasionally 
one observes systematic peak modulations as reported in 
Ref. ||. In striking contrast, the results for wide leads 
show a few isolated large peaks, separated by 15 — 25 
levels with much smaller peak height (not visible on the 
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FIG. 1. Coupling strength g p for a sequence of 100 states 
in the interval 290 < p < 390. Results are for a quantum 
dot with quadrupolar shape and deformation e = 0.2. (a) 
Narrow leads: ka e g = 0.1. (b) Wide leads: ka e g = 5.0. Here 
k = (2m*E p ) 1/2 /h is the wave number associated with the 
bouncing-ball state p — 337. Inset: Real-space projection of 
the state p = 337. 



scale of Fig. 1(b)). All large peaks are associated with 
states quantized on stable bouncing-ball orbits. This is 
illustrated in the inset for the state p — 337. The height 
of the small peaks not resolved in Fig. 1(b) is typically 
two or more orders of magnitude smaller than the maxi- 
mum peak height. Such tiny peaks are difficult to resolve 
in standard Coulomb blockade experiments. The inter- 
ference experiments face additional noise from the 
transmission through the free arm of the ring. For wide 
leads, interference experiments are therefore only sensi- 
tive to the strongly coupled bouncing-ball modes. 

The crucial role of the transverse barrier width a e g 
can be understood from the following argument: Con- 
finement in a lead of width a c g yields the transverse mo- 
mentum spread h/a e g for electrons injected in the quan- 
tum dot. Wide leads therefore result in near normal in- 
jection and provide exceptionally strong coupling to the 
bouncing-ball states. To quantify this argument we ex- 
press the partial width T l p — \V p \ 2 in terms of the Husimi 
function H p . Recalling the Gaussian form of ipi and us- 
ing the arc length s as integration variable in Eq. (^), we 
find 



r 



\ Cl w 



H p (s l ,0), 



(4) 



where H p is calculated for the normal derivative of the 
wave function ip p . The Husimi function is evaluated at 
the phase space point (s = s l ,p s = 0) reflecting the po- 
sition si of the lead and zero transverse momentum. We 
note that the Husimi function can be written as an inte- 
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gral over the Wigner function, smoothed with Gaussians 
both in arc length and transverse momentum. Their 
width is given by a c g/y/2 and h/y/2a c s, respectively. We 
express these widths in terms of the uncertainties Aip, 
Asinx in polar angle <f> and injection angle Xi respec- 
tively (x is the angle between the momentum of the in- 
jected electrons and the normal to the boundary). Using 
the relation p s = psinx, we find A(f> w 0.7(fca e ff)/(fci?o) 
and A sinx 0.7/(fea e g), where i?o is the average radius 
of the dot. 

The phase space portrait of the dot with deformation 
e = 0.2 is shown in Fig. 2 using (</), sinx) coordinates. 
Superimposed is the Husimi function of the bouncing- 
ball mode p = 337 calculated for ka e g = 5.0 and fci?o ~ 
40, corresponding to approximately 400 electrons on the 
dot. The Husimi function attains maximum value at the 
phase space points (</>, sinx) = ±(7r/2,0) representing 
normal injection from the leads. 

In order to obtain the Husimi function and hence T p 
analytically, we must solve for the wave function ip p . For 
states localized on a chain of stable islands the calcu- 
lation can be carried out using the semiclassical theory 
developed in Refs. |ll|[l4|]. The result can be expressed 
in terms of the stable periodic orbit at the center of the 
island, 



H p (s,p s ) = 



A„ 



As„Ap, 



■ exp 



x exp 



Ap, t 



(5) 



where s M is the arc length and p^ the transverse mo- 
mentum at the bounce points /j, of the periodic orbit. 
Moreover, As^ = [a 2 ff + Z 2 ] 1 / 2 a nd A Pfl = [(h/a cS ) 2 + 
{h/l^f] 1 ' 2 where l„ = y/2h\m% 2 \\4- Tr 2 ^!^- 1 / 4 is de- 



termined by the monodromy matrix M M = (m^). The 
amplitude depends only weakly on a e g . The result (|^) 
allows us to express the width of island states in terms 
of the phase space distance of the underlying periodic 
orbit from the lead injection points. No comparable ana- 
lytical result is known for states quantized in the chaotic 
sea. However, we find numerically that for wide leads the 
Husimi function of such states has exponentially small 
support near the injection points (see the inset of Fig. 2). 

We now turn to the calculation of the transmission co- 
efficient Tqd and the phase <^qd- We assume wide leads 
and temperature kT < A. The calculation for the weak 
coupling limit T p <C kT, A is straightforward and can be 
done using the Green function derived in Ref. jn]]; the 
results will be presented elsewhere Jl?]]. Here, we address 
the regime of a more open dot as realized in the experi- 
ment [|| . This regime is characterized by Tp ~ A where 
we identify p with the bouncing-ball state closest to the 
Fermi energy. All other states p ^ p in the vicinity of Ep 
have a width r p <C A. The Green function is found us- 
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FIG. 2. Husimi projection (color scale: red indicates high 
probability) of the bouncing-ball mode p = 337 calculated for 
ka g = 5.0 and e = 0.2. Superimposed is the Poincare surface 
of section of the classical phase space in {<f>, sin x) coordinates. 
Inset: H p (<j> — tt/2,0) vs. fea e ff for the states p = 337 (solid 
line) and p = 347 (dashed line, enhanced by a factor 18), the 
latter state is localized in the chaotic region of phase space. 



ing the equations-of-motion method [Q . It is diagonal 
up to small off-diagonal corrections 0(y / T p Tp/ A), and 
given by 



Gp p (E) - ^ ~£TZ 



Pn 



N=0 



E-[Ep- eV g + UN] + iTp/2 ' 



(6) 



Here N counts the total number of electrons in all dot 
levels except for the level p, Pn is the respective occupa- 
tion probability, Ep is the single-particle energy of state 
p, V g is the gate voltage, and U = e 2 /C is the charg- 
ing energy. The transmission through the states p p 
is negligible. Equation (Q) is the generalization of the 
weak-coupling result jllj to the case of a single strongly 
conducting quantum state. The probability P^ = 
Z- l exv[-Vt(N)/kT] with Z = £ N exp[~n(N) / 'kT] is 
related to the thermodynamic potential fl(N) of the dot. 
To evaluate Pn we replace ft n by fl N + [E p — eV g + UN — 
EF]{np) n, where Q N is calculated for the dot with level 
p excluded from the spectrum, and 

(np) N =-- fdElm = tlB. (7) 

V Pl nj E-[E p -eV g + UN]+iT p /2 w 

is the canonical occupation probability of level p. 

In Fig. 3 we show the transmission coefficient Tqd 
and the phase <Pqd vs. gate voltage V g calculated for 
kT = 0.2A, T p = 1.5A and U = 12A. All peaks shown 
result from transmission through the level p — 337. The 
peaks have comparable height and similar phase in quali- 
tative agreement with the experiments |ra,pl . The central 
peak has a Lorentzian shape of width Tp. Note that 
the transmission peaks develop a peculiar asymmetry as 
the conducting level moves away from the Fermi energy: 
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FIG. 3. (a) Transmission coefficient Tqd and (b) phase 
<^>qd as a function of gate voltage V g evaluated for kT = 0.2A, 
Tp = 1.5A, and U — 12A. The peak asymmetry in (a) is a 
unique fingerprint of BBT. 



Each peak to the left and to the right of the central peak 
falls off more rapidly on the side facing the central peak. 
This pattern extends over the whole sequence and be- 
comes more pronounced for the peaks in the tails. This 
asymmetry results from the breaking of the particle-hole 
symmetry in the transmission through the bouncing-ball 
state and is a unique fingerprint of BBT. Inspection of 
the data of the experiment H reveals the same asymme- 
try [jl9| , providing strong evidence that the peak correla- 
tions in this experiment are due to BBT. We note that the 
number of correlated peaks in Fig. 3 is less than observed 
in experiment. The origin for this is addressed below. A 
detailed discussion of this issue and of the phase lapse 
between the peaks will be given in a future publication 

0- 

We finally address the universality of the results pre- 
sented above. The stability of bouncing-ball motion does 
not rely on the effective dot potential considered here 
Similar stability if found for dots of other shapes [|l3 14 
as well as for smooth confining potentials ||. We as- 
sumed that the dot shape does not change upon the ad- 
dition of electrons. Some robustness of the shape has in- 
deed been demonstrated in a recent experiment ]2C| ], In 
general, however, variation of the gate voltage is expected 
P,pT| to change the electrostatic potential and modify 
the shape of the dot. It was shown |9|]21j] that shape de- 
formations can enhance peak correlations by "pinning" 
conducting levels close to the Fermi energy. To study 
the effect of deformations, we varied the shape by a pa- 
rameter linear in gate voltage V g . Over a range of V g 
we observed the predicted enhancement of correlations 



and sequences of more than 10 correlated peaks. At the 
same time, the bouncing-ball modes changed very little 
with V g as substantial modification of the shortest stable 
modes typically requires a large change in potential. 

In conclusion, we have demonstrated a new mecha- 
nism for transport through quantum dots in the Coulomb 
blockade regime. Current flows by tunneling through 
bouncing-ball modes. The fingerprints of bouncing-ball 
tunneling are sequences of correlated conductance peaks 
with asymmetric peak line shape. We find the peak 
asymmetries in recent Coulomb blockade interference ex- 
periments, providing strong evidence that the peak corre- 
lations in these experiments are caused by bouncing-ball 
tunneling. 

We acknowledge helpful conversations with F. Haake 
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the Deutsche Forschungsgemeinschaft, DGAPA-UNAM, 
and CONACyT-Mexico. 
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